Scarring by homoclinic and heteroclinic orbits 
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In addition to the well known scarring effect of periodic orbits, we show here that homoclinic and 
heteroclinic orbits, which are cornerstones in the theory of classical chaos, also scar eigenfunctions of 
classically chaotic systems when associated closed circuits in phase space are properly quantized, thus 
introducing strong quantum correlations. The corresponding quantization rules are also established. 
This opens the door for developing computationally tractable methods to calculate eigenstates of 
chaotic systems. 

PACS numbers: 05.45.Mt, 03.65.Sq 
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Scarring constitutes an important topic in the field of 
quantum chaos 0, Q. S S 0; IE ' representing a dra- 
matic departure from random matrix theory predictions 
0, and providing many practical applications j^. The 
term "scar" was introduced by Heller to describe an en- 
hanced probability density in certain eigenfunctions of 
classically chaotic systems taking place along short pe- 
riodic orbits (POs) P, I3- In order to explain this phe- 
nomenon, he noted that recurrences along a PO can com- 
pete under suitable conditions with the unstable dynam- 
ics in its vicinity, giving rise to the building up of prob- 
ability. From this point of view, an interesting specula- 
tion jj, 0] concerns the possibility that long time mo- 
tions, such as homoclinic or heteroclinic orbits, which are 
cornerstones in the classical theory of chaos, could also 
generate scars. 

In this Letter, we demonstrate that such motions can 
indeed produce enhanced probability density accumula- 
tion associated with POs. However, this new type of 
localization differs in two aspects from that previously 
described in the literature [jJ. In the case of scars pro- 
duced by POs, quantized closed circuits are defined in 
the direction of the dynamical flux, and attention is fo- 
cussed into the part of the probability remaining near 
the PO after recurrences. In our case, the relevant cir- 
cuits to be quantized are defined in a transverse direction 
to the flux, clearly observed on a surface of section, and 
the attention is concentrated on the probability leaving 
and then returning to the vicinity of the PO. This is a 
nontrivial fact, since it implies that the time scales char- 
acterizing these processes are of very different nature. 
While PO scarring takes place in times of the order of 
the PO period, the localization induced by homoclinic 
and heteroclinic excursions needs a time of the order of 
the Ehrenfest time, t^, to develop. 

The model used in our calculations is the fully chaotic 
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FIG. 1: (a) Tube and (b) scar wave functions along the hori- 
zontal periodic orbit of a desymmetrized stadium billiard cor- 
responding to the quantum number uh = 44. The associated 
Husimi based quantum surfaces of section, in Birkhoff co- 
ordinates, are shown in (c) and (d). The manifolds of the 
horizontal periodic orbit are also plotted in panel (d). 



system consisting of a particle confined in a desym- 
metrized Bunimovich stadium billiard, defined by the ra- 
dius of the circular part, r = 1, and the enclosed area, 
I+tt/A. In the quantum calculations, Dirichlet conditions 
on the stadium boundary and Neumman conditions on 
the symmetry axes are imposed. We focus our atten- 
tion on the dynamics influenced by the horizontal PO. 
To study scarring phenomena related to this orbit, we 
deflne scar functions in the following way. We start from 
"tube functions", |(/'tubo), covering the region around the 
PO, obtained with the semiclassical theory of resonances 
developed in Ref. 0. Associated to the construction of 
these wave functions there is a quantization on the PO 
energy, i?BSi given by the usual Bohr-Sommerfeld (BS) 
condition on the wavenumber 
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being uh the quantum number along the orbit, Lh = 4 
its length, and uh ~ 3 the corresponding Maslov in- 
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FIG. 2: Husimi based quantum surfaces of section for the 
scar functions corresponding to the quantization condition 
summarized in Table |l] showing the localization effects pro- 
duced by the homoclinic motion on hoi and ho2, when the 
corresponding areas (see Fig. |^ are quantized. To guide the 
eye we have superimposed the manifolds corresponding to the 
horizontal periodic orbit. 



dex. This condition guarantees the coherent contribu- 
tion of probability recurrences along the PO. An exam- 
ple of these tube functions is shown in Fig.^a). When 
examined in phase space, using a Husimi based quan- 
tum surface of section (QSOS) Il3|, this function appears 
localized around the associated classical fixed point at 
{Q:P){A) — (1 +7r/2, 0); see Fig. ^c). In a second step 
the dynamical information up to a given time, T, is in- 
corporated into these functions by using the dynamical 
averaging method developed in Ref. fl^ 
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where the cosine function has been included in order to 
minimize the energy dispersion of |</>scar) ^^^^ 
way, all processes transporting probability from and to 
the region defined by \4>tuhc), in times smaller than T 
and with energies in the window E-Qg ± h/T, will be 
taken into account ^Jl- An example for T — 0.95tE is 
shown in Figs.^^b) and (d). As can be seen, a self-focal 
point appears in the wave function [Fig. H^b)] and the 
phase space density spreads along the hyperbolic struc- 
ture as the manifolds go away from the PO [Fig. H^d)]. 
We emphasize that in this case, the only channel able to 
transport probability from and to the region defined by 
I'/'tubc) is the PO, this being the reason that a semiclas- 
sical recipe for the construction of scar functions can be 
easily derived 

Let us now consider longer propagation times, thus 
opening the possibility to observe more interesting dy- 
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40 


64.010 


34.07 30.47 


(c) 


44 


70.293 


37.43 


33.46 


(d) 


50 


79.718 


42.49 


37.95 



TABLE I: Quantum number and Bohr-Sommerfeld quantized 
wavelength along the horizontal PO [Eq. , and homoclinic 
quantum numbers [Eq. (|2J] corresponding to the cases shown 
in Fig. H 



namical effects. In Fig. |21we show QSOS for some scar 
functions calculated by means of Eq. Q with T = lAOts, 
so that the effects of the primary homoclinic motion asso- 
ciated to the horizontal PO have had time to fully develop 
(notice the different range spanned by the vertical axes 
with respect to Fig. . The different panels correspond 
to four BS quantization conditions on the horizontal PO, 
which are summarized in Table ^ As can be seen, all 
QSOS show a prominent peak on the horizontal PO fixed 
point, with much of the probability density spreading all 
along the full length (linear and nonlinear parts) of the 
emanating manifolds. Furthermore, a closer examination 
reveals conspicuous accumulations of density on different 
points along these manifolds. To analyze in detail which 
are these points, we present in Fig. 13 a detailed phase 




FIG. 3: Detail of the phase space portrait for the desym- 
metrized stadium billiard in Birkhoff coordinates. 
(Left) Fixed point, (A), and unstable (full line) and stable 
(dashed line) manifolds corresponding to the horizontal peri- 
odic orbit. These manifolds first cross, with different topol- 
ogy, at the primary homoclinic points, hoi and ho2, which 
map (backwards in time) into ho'i, . . ., and ho2, . . ., respec- 
tively. They define two primary homoclinic areas, Shoi and 
Sho2, which are shown in the inset. 

(Center) Fixed point, (B), and manifolds for the "flipped L"- 
shaped periodic orbit (see inset). The heteroclinic points with 
the horizontal orbit, he'i and he* , and the corresponding het- 
eroclinic also shown. 
(Right) Same for the "triangle" orbit. 
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space portrait of the classical structures related to the 
scarring orbit. In the left part we see the stable and un- 
stable manifolds emanating from fixed point (A), which 
first cross (with different topology) at the two primary 
homoclinic points hoi and ho2- These points map (back- 
wards in time) into the sequences: ho'^, . . . and /lOj, . . ., 
respectively, which constitute the primary homoclinic or- 
bits defining two relevant homoclinic areas, Shoi and 
Sho2j in phase space. Now, comparing this with Fig.[21[a) 
we see that the density maxima in the QSOS appear lo- 
cated at the two primary homoclinic points hoi and ho2- 
If we continue examining Fig.|2 it is observed in (b) that 
only the maximum on hoi exists, presenting the QSOS 
a minimum of density on the other primary homoclinic 
point, ho2- For case (d) the situation is reversed, appear- 
ing a minimum on hoi and a maximum on ho2- Finally, 
we find in (c) a case with two minima at hoi and ho2, 
respectively. 

This behavior can be explained by taking into account 
that, associated to the homoclinic orbits, there are fluxes 
of quantum probability, which follow the corresponding 
circuits in phase space. Each of these homoclinic excur- 
sions can be thought of as a channel for this flux. In 
these excursions, some part of the quantum probability 
leaves the vicinity of the PO and returns later to it with 
some phase, shifted with respect to the phase accumu- 
lated by the flux along the PO circuit. When this phase 
shift is a multiple of 27r, the two probability fluxes in- 
terfere constructively. In this case, the channel can be 
considered as a bright channel for the process, while in 
the case of a destructive interference the channel will be 
dark. The location of these channels is best ascertained 
by examination of the QSOS, where the density at the in- 
tersection with the homoclinic orbits is monitored. When 
one channel (homoclinic orbit) is bright an accumulation 
of density on the corresponding homoclinic points will be 
visible, while for the dark channels the density at those 
points is minimum. Furthermore, the phase shifted by 
an homoclinic excursion consists of two different contri- 
butions. One, of dynamical origin, is given by the area 
enclosed by the homoclinic circuit, while the other is of 
topological nature. Correspondingly, a constructive in- 
terference is governed by the following quantization rule 
derived in Ref . Il7| 



kSho = 27r(?i,io + t^/io/4). 



(3) 



where Uho 



1 and i^ho 



in the case we are con- 



sidering here. In this way, the accumulations or lack of 
probability density observed in the plots of Fig. [3 can be 
explained as the result of the quantization and/or anti- 
quantization of the primary homoclinic areas, Shoi or 
Sho2 ■ This is numerically illustrated in Table where 
the associated "homoclinic" quantum numbers, rihoi and 
niio2, computed at the quantized wavelength values of 
the horizontal PO, A:bSj (taken as approximations to the 
true khoi 2 values) are given. Moreover, the four cases 
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FIG. 4: (a) Husimi based quantum surface of section (QSOS) 
for the horizontal scar function with n — 224 and T — 3.3tE- 
The inset shows the intensities in the basis set of the stadium 
eigenfunctions. (b) Profile of the QSOS along the unstable 
manifold (full line). The results for T — OMe (dotted line) 
and T = 1.2tE (dashed line) are also included. 



selected correspond to situations in which Uhoi 2 ^^^y 
very good accuracy, either integer (quantization) or half- 
integer (anti-quantization) numbers. 

Notice, that in the previous discussion only the two 
primary homoclinic orbits, defining the shortest possi- 
ble homoclinic circuits, have been used. The question 
now arises as to what is the influence of other homoclinic 
orbits, performing longer and more complicated excur- 
sions into phase space regions where the dynamics are 
influenced by other POs jSj] . These secondary orbits also 
define phase space circuits that can be quantized in a way 
similar to |j2Jl. To investigate their role in the issue ad- 
dressed in this Letter, we show in Fig.^a) the QSOS for 
the horizontal scar function with jih — 224, computed 
for T = S.Me- In this case, both primary homoclinic 
areas are quantized, with values of rihoi — 189.01 and 
n/io2 = 168.07. As can be seen, a lot of fine structure has 
developed here and several density accumulations appear 
along the manifolds. To analyze more easily the result- 
ing structure, we have plotted in (b) the QSOS profile, 
H[q,pm{q)], along the unstable manifold of the horizon- 
tal PO. Three values of the evolution time have been 
considered, corresponding to T/ts = 0.9, 1.2, 3.3, re- 
spectively. For the shortest value (dotted line), only the 
maxima on the primary homoclinic points, hoi and /102, 
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are (barely) visible. This result is in agreement with our 
previous conclusion about the results of Fig. Ufa). As 
the evolution time increases, and we move to T = 1-21^ 
(dashed line), the structure is better resolved. As a con- 
sequence, the two previous maxima get more prominent, 
and also a new peak centered on he2 appears. Finally, 
at the longest time considered, T = 3.3tE, all previous 
maxima get better defined, two extra maxima appear on 
hei and he'i, and even successive mappings of previous 
points, namely ho'i, are resolved within the big peak at 
(Zm = 1 + 7r/2 [corresponding to (A) in Fig. Following 
our previous discussion, the existence of these new peaks 
could be explained as a consequence of the quantization 
of the phase space circuit defined by secondary homo- 
clinic orbits. Rather, we prefer to use here an argument 
based on heteroclinic orbits. The advantage being that in 
this way a family of homoclinic orbits leaving the original 
PO, making an integer number of turns around a new vis- 
ited short PO, and returning again to the vicinity of the 
original PO, can be substituted by a single heteroclinic 
circuit [ij, [l^ . This circuit consists of an heteroclinic 
orbit that populates the new PO by transporting proba- 
bility from the original PO, and another heteroclinic orbit 
transporting probability in the opposite direction. In this 
way, long term dynamics are dealt with more efficiently, 
thus avoiding the serious computational problem associ- 
ated to the exponential proliferation of homoclinic orbits. 
These heteroclinic circuits are identified, on a surface of 
section, by four consecutive pieces of manifold, defining 
an invariant (heteroclinic) area in phase space. For in- 
stance, the heteroclinic circuit between the horizontal PO 
and the "flipped L" -shaped (rendering an hexagon in the 
full version of the stadium) PO (see Fig.|3|center) is given 
by: the piece of unstable manifold from (A) to he[, the 
piece of stable manifold from he'i to the (B), the piece of 
unstable manifold from (B) to the heteroclinic point he* 
(belonging to the heteroclinic orbit transporting proba- 
bility towards the horizontal PO again), and finally, the 
piece of stable manifold from he* to (A) . The transverse 
area enclosed by the circuit (shaded in the figure) is Shei ■ 
We emphasize that by evolving this circuit forward and 
backward in time, an heteroclinic tube characterized by 
Shei is defined. Having this in mind, the point he[ and its 
mapping hei are assigned to the same heteroclinic chan- 
nel, which will be bright when the following quantization 
condition is fulfilled, 

kShei = ^imhei- (4) 

This condition, derived from the fact that the overlap 
between the scar function for T > 2tE and the fiipped-L 
tube function is proportional to cos{kShei/2) 14], guar- 
antees that the region around the flipped-L PO is ef- 
fectively populated through the heteroclinic channel de- 
scribed above. In our case, rihei — 19.00, which is ex- 
tremely close to an integer value, thus explaining the exis- 
tence of the corresponding maxima in the QSOS. Another 



maximum appear at he2 in Fig. 0] It can be explained 
through the heteroclinic circuit associated to the "tri- 
angle" ("kite" in the full stadium) PO and related area 
5/162 shown in the left part of Fig. 13 In this case the 
computed value for the "heteroclinic" quantum number 
is n/iea — 5.98, again extremely close to an integer value. 

A final point worth discussing concerns the spectral lo- 
calization of the scar functions ((SJ, that have been used 
throughout this Letter. For this purpose, an inset with 
the squared modulus of the coefficients projecting the 
corresponding scar function on the stadium eigenfunc- 
tions has been included in Fig. 0] As can be seen, it es- 
sentially consists of only 7 contributions, a terribly small 
number if one takes into account that a generic plane 
wave (usually taken as optimal basis functions for the 
evaluation of billiard eigenfunctions) with the same k re- 
quires at least 300 eigenfunctions for a correct descrip- 
tion. This clearly indicates that the scarring effect by 
homoclinic and heteroclinic orbits that we have unveiled 
here really opens the door for a systematic construction 
of the eigenfunctions of classically chaotic systems. 
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